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The Lya forest in quasar spectra probes scales where massive neutrinos can strongly suppress the 
growth of mass fluctuations. Using hydrodynamic simulations with massive neutrinos, we success- 
fully test techniques developed to measure the mass power spectrum from the forest. A recent 
observational measurement in conjunction with a conservative implementation of other cosmolog- 
ical constraints places upper limits on the neutrino mass: m„ < 5.5 eV for all values of Q m , and 
m v < 2.4(fi TO /0.17-l) eV, if 0.2 < fl m < 0.5 as currently observationally favored (both 95% C.L.). 



Experimental evidence for finite neutrino masses and 
flavor oscillations continues to mount. Recently the 
Super-Kamiokande experiment has provided strong ev- 
idence that oscillations from to another species in- 
volve a mass greater than V dm 2 = 0.07~^^eV jjj. The 
LSND experiment suggests the existence of to v e os- 
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0.003 eV Si. These 
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neutrino deficit requires V Sm 2 
mass splitting results are consistent with one to three 
weakly interacting neutrinos in the eV mass range Q . 

Neutrinos in this mass range are important cosmolog- 
ically since, if they exist, they would represent a non- 
negligible contribution to the dark matter content of the 
Universe. In units of the critical density, neutrinos con- 
tribute fl v — Nh~ 2 m u /%A eV, where h is the dimension- 
less Hubble constant (Ho — lOO/i km s _1 Mpc -1 ), and 
N is the number of degenerate mass neutrinos. As light 
neutrinos do not cluster on small scales, they retard the 
gravitational growth of density fluctuations. Any mea- 
sure of small-scale clustering is thus sensitive to neutrino 
masses in this range. One such useful measure is the 
clustering of the intergalactic medium revealed by the 
absorption features in quasar spectra known as the Lya 
forest j!| . As we will discuss below, the Lya forest has the 
distinct advantage that clustering properties of the mass 
distribution can be inferred from it, which greatly facili- 
tates comparisons with theory, and makes large searches 
of parameter-space possible. 

In this Letter, we make use of a recent Lya forest mea- 
surement of the power spectrum of mass fluctuations || 
to place limits on the mass of the neutrino(s). Conver- 
sion of the power spectrum measurements into neutrino 
mass limits requires a framework for cosmological struc- 
ture formation. There is growing evidence that struc- 
ture formed by the gravitational instability of cold dark 
matter (CDM) with adiabatic, Gaussian initial fluctua- 
tions. Upcoming Cosmic Microwave Background (CMB) 
experiments should conclusively determine whether this 
assumption is a good one (7|. For now, we note that 
this framework, which we adopt, includes all currently 



favored models, and also that whatever the true model 
for structure formation, Lya forest measurements should 
respond sensitively to the presence of massive neutrinos. 

Our adiabatic CDM dominated universes are described 
by 6 free parameters: the matter density f2 m , dimen- 
sionless Hubble constant h, baryon density Sib, neu- 
trino density f2„, density fluctuation amplitude A, and 
tilt n, which define the initial density power spectrum 
-Pinit(fc) = Ak n . We also initially assume that spatial 
geometry is flat, as implied by recent measurements of 
distant supernovae and CMB anisotropies (|], before in- 
vestigating the consequences of relaxing this assumption. 

We begin by describing the Lya forest power spectrum 
measurement method and test it on hydrodynamic sim- 
ulations. We then apply the observational constraint to 
the 6 dimensional CDM parameter space to find an upper 
limit on the neutrino mass. Given this large parameter 
space, we conservatively employ other cosmological con- 
straints, notably from the abundance of galaxy clusters 
and the age of globular clusters, to constrain other pa- 
rameters that can mimic the effects of massive neutrinos. 
Finally, we consider prospects for making a precise mea- 
surement of m„ using future Lya forest observations and 
upcoming CMB experiments. 

Testing Lya forest simulations with m v > 0. — The 
Lya forest of neutral hydrogen absorption seen in quasar 
spectra |5j arises naturally in cosmological scenarios 
where structure forms by the action of gravitational in- 
stability. In hydrodynamic simulations of such models |J 
most of the absorption arises in gas of moderate overden- 
sity, whose physical state is governed by simple processes 
(mainly photoionization heating and adiabiatic cooling, 
see also the analytical modeling of (l(J ) . The density field 
can then be locally related to the optical depth for Lya 
absorption |p"l|| and hence a directly observable quantity, 
the transmitted flux in a quasar spectrum. 

The Lya forest can therefore be used to determine the 
statistical properties of the density distribution, and in 
particular P(k), the power spectrum of density fluctua- 
tions. A method for carrying this out was described by 
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|L2| , who also tested it on hydrodynamic simulations with 
CDM only. The method relies explicitly on the assump- 
tions that the initial fluctuations were Gaussian, and that 
gravitational instability was responsible for their growth. 
A measurement from an observational dataset was made 
by As we will use this result to constrain the neu- 
trino mass, we first test the method on a hydrodynamic 
simulation which includes massive neutrinos. 

The measurement of P(k) from Lya forest spectra is 
carried out in two stages. First, the shape of P(k) is mea- 
sured from the power spectrum of the Lya forest flux. 
Second, normalizing simulations are used to set the am- 
plitude of the linear mass P(k). We refer the reader to 
@ for details. 

The hydrodynamic simulation itself is described in de- 
tail by [|13|. We follow the evolution of structure in a 
model with two mass-degenerate neutrino species, us- 
ing the Parallel TreeSPH hydrodynamic code jwj. The 
model parameters are Q m = 1, h = 0.5, = 0.075, and 
f2„ = 0.2, so that the mass in both species combined 
is 5 eV. This so-called cold plus hot dark matter model 
(CHDM) is normalized to fit the COBE results p), so 
that the amplitude of mass fluctuations in 8 h~ 1 Mpc 
spheres at z = 0, a$ — 0.7. We use a box of size 
11.111 /i Mpc, periodic boundary conditions and ini- 
tial conditions taken from |l6(| . The CDM and gas com- 
ponents are represented by 64 3 particles each, and the 
neutrinos by 2 x 64 3 particles. We use the distribution 
and physical state of the gas at z = 2.5 to generate artifi- 
cial Lya spectra, for 1200 randomly chosen lines of sight 
through the simulation volume. 

We then apply the P(k) recovery method of [Q to 
these spectra. We use normalizing simulations run under 
the PM approximation 0,^2| with 64 3 particles and an 
ll.lllft,^ 1 Mpc box. We use the same estimator for the 
amplitude of P(k) as in the observational analysis paper, 
|[| . The results of the test are shown in Fig. 1, where we 
plot the recovered P(k), together with the linear theory 
prediction for the model. We also show the linear power 
spectrum of a CDM-only model (with O m = 1, h = 0.5), 
again normalized to COBE, so that a% = 1.2. 

The error bars on the P(k) points are representative of 
the "cosmic variance" error which arises from only having 
one hydrodynamic simulation volume. We estimate this 
uncertainty by running 10 additional PM approximation 
simulations of the CHDM model, and extracting 1200 
lines of sight from each. The standard deviation of their 
results for each P{k) point provides the error bar. There 
is an additional overall amplitude uncertainty associated 
with the normalization. We estimate this by applying 
the normalizing procedure to the PM CHDM simulations 
taken individually, finding that the additional error in 
this test case is a negligible 4% in P(k). 

We use the simulation to test for systematic errors in 
the technique. The observational result was given by || 
in terms of a power law fit to the P(k) data points with 
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FIG. 1. A test of the P(k) recovery method. The points 
are P(k) recovered from a hydrodynamic simulation of the 
Lya forest in a massive neutrino model (CHDM see text). 
The turndown for k > 5/i.Mpc -1 is caused mainly by the lim- 
ited resolution of the simulation and lies outside of the obser- 
vational band (shaded). The lines show the linear P(k) pre- 
dictions for the massive neutrino model, and for a CDM-only 
model. 



2.7 x 10 3 (kms" 1 ) < k < 1.42 x 10 2 (kms" 1 ) (which 
corresponds to 0.5 hMpc^ 1 < k < 2.7 /iMpc" 1 for an 
Einstein-de Sitter model). We fit the simulation data 
points to a power-law over this range, finding an ampli- 
tude A 2 (ft) = 0.16 ± 0.028 (la) at ft = 1.5/iMpc" 1 [= 
0.008 (kms" 1 ) -1 ]. Here A 2 (k) = k 3 P(k)/2n 2 , the contri- 
bution to the density field variance from a unit interval 
in log ft. The logarithmic slope, n = — 2.18to'28 (l* 7 )- 
The linear theory prediction for CHDM is A 2 (ft) = 0.21, 
n = —2.40. Our recovered P(k) is therefore about 2a too 
low in amplitude, and has a slightly flatter slope. By ex- 
amining results from the more numerous PM simulations, 
we find that the largest scale data point is systematically 
lowered by peculiar velocity distortions (as predicted by 
p8[ ) , an effect which is not accounted for in our estimate 
of the P(k) shape. Including or leaving out this point 
(which has the largest statistical errors) has only a small 
effect on the power law fit. It is possible, however, that 
taking these effects into account or further refining the 
analysis would improve the result. For the moment it is 
sufficient to note that the observational result of || cur- 
rently has a statistical uncertainty of ~ 70% (2a), larger 
than any biases revealed by our test. 

Constraining the neutrino mass. — Although the sup- 
pression of power in the Lya forest due to a finite neu- 
trino mass is large (~ 800ft,y/n m percent [ft9|), other as- 
pects of the cosmological model can counterbalance this 
effect; limits on the neutrino mass consequently depend 
on the range of models allowed by other cosmological 
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constraints. 

In addition to the Lya forest power spectrum measure- 
ment of A 2 (0.008[ km s" 1 ]- 1 ) = 0.57l{j;l$ (95% CL) (we 
do not use the Lya slope measurement, which has no 
significant effect) at z = 2.5, we consider 6 other cosmo- 
logical constraints. The Hubble constant is measured to 
be ft. = 0.72±0.17 (95% CL), where we have added statis- 
tical and systematic errors in quadrature and doubled the 
lcr errors p0| |. The amplitude of the fluctuations is deter- 
mined by the COBE detection of large angle anisotropics 
|L5[ ; we ignore the 7% measurement uncertainties on the 
temperature fluctuations which are substantially smaller 
than the other uncertainties. The abundance of galaxy 
clusters today constrains models at the 8/i _1 Mpc scale, 
where the amplitude of fluctuations in top-hat spheres 
(at 2=0) is a 8 = 0.56^ 47 with 95% CL of +20^ 21og n - 
percent and -18ft^ 21ognm percent |2l). 

Galaxy surveys measure the shape of the power spec- 
trum; we use the compilation of |p2|| and take measure- 
ments only from the range 0.025/i Mpc -1 < k < 0.25/i 
Mpc -1 to avoid spurious survey volume effects on the 
large scale and uncertainties in the nonlinear corrections 
on the small scale. For the galaxy survey data, we employ 
a Ax 2 statistic and take A\ 2 = 4 to represent the 95% 
confidence limits on the shape of the power spectrum. In 
carrying this out, we assume that a linear "bias" is oper- 
ating so that the matter power spectrum is related to the 
galaxy power spectrum by a constant multiplicative fac- 
tor (which we find as part of our \ 2 minimization). We 
also employ nucleosynthesis constraints on the baryon 
density of fl b h 2 = 0.019 ±0.0024 (95% CL) @. Finally, 
we place a lower limit on the age of the universe by as- 
suming that it must at least as old as the oldest globular 
clusters (13.2 ± 2.9)Gyrs (95% CL) @). 

Given these constraints, we could construct a joint like- 
lihood to find the best fitting neutrino mass. We have 
decided to be more conservative, however, and consider a 
model ruled out if it violates the (95%) confidence limits 
on any constraint taken individually. With these "2c" 
constraints on the parameter space, we use the analytic 
approximations of [ ^6| to explore the remaining space 
rapidly and find the model that maximizes the neutrino 
mass as a function of Q m . The result for one massive 
species is displayed in Fig. ||a. We also show the effect 
of omitting the Lya forest measurement. The measure- 
ment has a powerful constraining effect at low and high 
Q m since the amount of tilt required to match the clus- 
ter abundance and galaxy power spectrum shape violates 
the upper and lower Lya forest bounds respectively. For 
no value of Q m , can m„ be greater than 5.5eV. 

To address the robustness of our upper limits, we show 
the effect of scaling all errors by a factor of 1.5 to ap- 
proximate "3<t" constraints in Fig. 0. We also test for 
single point failures by dropping each constraint sequen- 
tially. Omission of either the age or cluster abundance 
constraint changes the maximal neutrino mass to ~ 7eV. 
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FIG. 2. Constraints on the neutrino mass (a) single massive 
species (b) two (degenerate) mass species with and without 
the Lya constraint. 



While dropping the galaxy power spectrum constraint 
does not increase the maximal mass substantially, it does 
weaken the bounds by up to 2 eV for f2 m < 0.5. Omission 
of the h or VL^h 2 constraint has essentially no effect. 

In applying our constraints, we assume that the uni- 
verse is flat (fi m + Qa = 1) and that gravity waves do 
not contribute to the COBE normalization. However, as- 
suming an open universe or gravity waves from power law 
inflation does not change the m v limits significantly since 
the tilt can be used to offset small changes in normaliza- 
tion. The simplest inflationary models can also predict 
a variation of the spectral index with scale. This is a 
small effect compared to the neutrino power suppression, 
so that we do not include it. Future CMB observations 
should address this point definitively. 

We also show the results assuming 2 neutrino species 
with identical masses in Fig. ||b. These limits are roughly 
half the single species results since the change in the 
growth rate is mainly governed by f2„. However dividing 
the total mass into 2 species makes each species more 
relativistic and enhances the suppression of the power on 
scales relevant to the galaxy power spectrum and cluster 
abundance constraints (27j. 

Future prospects. — Additional Lya forest data from 
quasar surveys such as the Sloan Digital Sky Survey 
(SDSS [p5[) have the potential to increase the precision 
of these mass constraints substantially. Furthermore, we 
expect that the next generation of CMB satellites will not 
only verify the existence of the underlying framework for 
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structure formation, which we currently assume, but also 
provide limits on the neutrino mass itself. 

How precise must these Lya measurements be to im- 
prove on projected CMB limits on the neutrino mass? 
To answer this question, we employ Fisher informa- 
tion matrix techniques to approximate the joint co- 
variance matrix. For the Lya forest power spectrum 
measurement, the Fisher matrix is given by F.y = 
(AP/P)- 2 (dlnP/dp i )(dlnP/dp j ). We add this to 
the CMB Fisher matrix projected for the MAP and 
Planck satellites including polarization information p8[ . 
The variance of the optimal unbiased estimator of pi 
marginalized over the other parameters is (F 

A fractional error of AP/P = 0.1 (la) on the Lya 
power spectrum would improve the MAP upper limit 
from 1.1 eV to 0.54 eV and the Planck limit from 0.51 eV 
to 0.29 eV both at 2a. Note that these represent limits 
in a wider 10 parameter space including spatial curva- 
ture and gravity waves. These improvements would ex- 
ceed those that can be achieved the SDSS galaxy sur- 
vey [p8| . CMB polarization information here is criti- 
cal for these improvements since they rely on an abso- 
lute normalization of the power spectrum from degree 
scale anisotropies. Polarization information eliminates 
the degeneracy between the normalization and the opti- 
cal depth due to reionization. 

Is 10% precision in power achievable from Lya for- 
est measurements? The constraint in this paper relies 
on data from ~ 10 full quasar spectra ||. The SDSS 
quasar survey J25| ] will yield spectra of roughly similar 
quality (resolution 2.5 A compared to ~ 1.5 A for 0]) 
for ~ 10 5 quasars. The mean distance between sight- 
lines in the SDSS will be somewhat larger than the scale 
on which clustering was measured by [0], so that the 
decrease in statistical errors will not be too far off the 
factor ~ (10 5 /10) 1 / 2 ~ 100 which would be expected if 
the sightlines were independent. A la statistical error of 
AP/P < 1% should therefore be possible in the future, 
so that systematic errors will become dominant. Study- 
ing larger hydrodynamic simulations should enable us to 
understand these systematic effects and refine our anal- 
ysis techniques. The vast size of the SDSS dataset will 
also enable us to pin down non-gravitational contribu- 
tions to Lya forest clustering, for example by analyzing 
the evolution of the forest with redshift. The signature 
of massive neutrinos, if they are present, should therefore 
be obvious, even if m v is a fraction of an eV. 
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